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Abstract 

Using chains of bistable springs, a model is derived to investigate the 
plastic behavior of carbon nanotube arrays with damage. We study the 
preconditioning effect due to the loading history by computing analytically 
the stress-strain pattern corresponding to a fatigue-type damage of the 
structure. We identify the convergence of the discrete response to the 
limiting case of infinitely many springs, both analytically in the framework 
of Gamma-convergence, as well as numerically. 
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1 Introduction 

Because of their interesting combination of properties, including high strength, 
low density, and high electrical and thermal conductivities, carbon nanotubes 
(CNTs) have been of great interest as nanoscale elements in a variety of applica- 
tions, [1]. Aligned arrays of CNTs can be readily synthesized to form foam- like 
materials that combine low density with a desirable dissipative response, [6l [10] . 
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There are a few interesting structural and mechanical features of these ma- 
terials: First, the thermal chemical vapor deposition process that is typically 
used to synthesize the arrays [B] results in a gradient in physical properties (such 
as density) along the height of the structure. This leads to strain localization 
during compression, with the majority of the structure remaining undeformed 
while increasing strain results in the sequential addition of highly localized buck- 
les [HI [11]. Second, there is a large amount of strain recovery after compression 
(typical of the CNT arrays that we study, which are synthesized using a vapor 
phase catalyst, [B], but not for all types of CNT arrays that are synthesized dif- 
ferently, [TH1[3]). Third, the stress-strain response of the material is hysteretic, 
with different loading and unloading paths. It has been noted in the past that 
repeated compressive cycles result in a hysteresis of decreasing area and decreas- 
ing stress at any given strain for the first few cycles, [6]. This effect, sometimes 
referred to as preconditioning, ceases after the initial few cycles, resulting in 
a hysteresis of constant area for loading cycles thereafter. Finally, we observe 
that this preconditioning effect is dependent on the maximum strain reached. 
When the maximum strain of all previous loading cycles is exceeded, the stress 
response obtained at these elevated strains is that of the un-preconditioned ma- 
terial, as if it had never been compressed previously, [O] . 

The features we have enumerated for the foam-like response of CNT arrays 
are analogous to those observed in other soft materials such as filled rubbers. 
In the context of rubbers, this response is associated with what is termed the 
Mullins effect, [T3j [8] . Capturing these features simultaneously in models has 
proven difficult in the past for other materials, with frequent use of simplifying 
assumptions that only allow models to match some of the experimental obser- 
vations (see, e.g., the idealized Mullins effect modeled in [5). 

We show in the present work that a suitable generalization of the mesoscopic 
mass-spring model of CNT structures recently proposed by ^ is able to handle 
the material damage due to preconditioning and permanent deformation within 
an effective one-dimensional framework. We introduce preconditioning-induced 
material damage by setting to zero the stiffness of a suitable percentage of the 
bistable springs that describe the response of the material at the microscopic 
scale. This allows us to model permanent axial deformation of the structure 
through the irreversible 'annihilation' of the springs with zero stiffness. 

2 Multiscale mass-spring models of CNT arrays 

2.1 Bistable spring model at the microscopic scale 

We model an infinitesimal portion of a CNT foam through the bistable spring 
model described in [9], which we hereafter briefly summarize. We assume that 
such a portion of the foam can be described as a chain of A'^ -I- 1 lumped masses 
mP, . . . , m^, with mP clamped at the bottom of the chain. The adjacent masses 
are connected to each other through bistable springs characterized by the axial 
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where u]y the axial displacement of the mass to' (positive upward), Hn 



L/N is the equal spacing between the masses, and un '■= {u%, . . . , 
potential and stress ct' vs strain laws of the generic spring are 



}. The 




:C2-fc^[£*-4+ln(l-(£'-4))], 
jLi e' i ^ i 

al + kl{e^-el), e\ < < t,, 

l-(e'-el)' 



£ < £ 

< e« < s' 

£^ < £% 



(2) 



(3) 



£' < £' 



where fcg > 0, fc^ < 0, A:* > 0, eJ^ > and £* > £„ are constitutive parameters 
(five independent parameters); the constants ci < and C2 > are such that 
V^iel,) = V^'isi), Vfc*(rj = V^isi); and it results (compare with Fig. 2 and Fig. 3 
of [9] for the notation) 



(4) 



with = /c^^, al^al + fc^(£l - £^). 
2.2 Plasticity and damage 

At the microscopic scale, the bistable springs introduced above permit a dy- 
namic switching process between the phases (a) and (c) , cf. [TH [TS] . As in [T3] , 
we name a response of the material plastic, if the strain £* of a single spring 
exceeds the threshold £^. For a chain of N springs, this can be characterized by 
the occurrence of loading and unloading stress plateaux. 

Within the current section, we rescale for simplicity L to unity; name (b) 
the unstable phase; and regard a mesoscopic element of a CNT array as the 
limit — > oo of a series of A^ microscopic springs. 

Let TO denote the number of hysteresis cycles that have been applied to the 
material during the previous loading history (up to different maximum strains). 
We assume that such a loading path has severely weakened the stiffness of 
(1 — P{m))N microscopic springs, for given < /3(m) < 1, and that at the 
current time it holds for all to £ N 

(Al) kl = kl for alH e N 

(y42) fcp = /cq = . . . = /cq = fco, fcg ^ — ... — k^ — 5 

{Ai) £:=£, for i = l,...,/3Ar 

(A4) £;=£j, for i = /37V + l,...,7V. 
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Condition (Al) stipulates the symmetry of the microscopic springs. (A2) states 
that the springs [pN + IJ to have stiffness fcg = 3, where (5 > is a smaU 
constant (damaged springs). We name undamaged those springs with stiffness 
constant fcg (springs 1 to PN). For the present analysis we also require a certain 
smallness condition on and £* relating to strong pinning that disappears for 
N ^ oo, see [Ti] , 

With given by ([2]), the mechanical energy of the structure is 



1 ^ 



Let (T be the given total stress. The mesoscopic average strain is simply e{uN) ■= 
^j^-^ e'(uAr), where e* denotes the strain associated with the i-th spring. 
Following an original idea of [15] , we model plasticity by the gradient flow 



equations 



iye\uN) = -^^(e^(uAr), . ■ .,e^{uN)) (5) 



with the total energy 

1 ^ 



The evolution equation ([S]) lets e* evolve towards local minimizers of $jv- We 
are interested in the limit f ^ which amounts to infinitely fast evolution such 
that e{uN) attains a local minimizer of $Ar. First we construct the equilibrium 
points. Inside the i-th spring element, the strain must satisfy the condition 

For given total stress a, there are at most the three local minimizers (using 
(A3)) 

£l{rn) = TT^, elim) ^ ^—p- + 4, (6a) 

Kg + (7 

^ ^ TT^ — = ^a("i) + £*• (6b) 



Note that for the derivation of we require that 5 is positive. 

In a loading or unloading experiment, the first spring located closer to the 
bottom of the structure is the softest and yields first, changing its phase, [16]. 
Next, the second spring yields, and so forth, until the /3A^-th spring. (Note that 
in accordance with (A4), the springs /3A^ + 1, . . . ,N with small spring constant 6 
do not flip.) Therefore, similar to the case of N identical springs, the total state 
of the series of springs is still completely specified by two scalar parameters p 
and q and the additional parameter (3. Here, p, q, 1 — p ~ q denote the phase 
fractions of the minimizers a, 5, and c, which corresponds to having jSNp, /3Nq 
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and (3N{1 — p — q) springs in phase a, 6, and c, respectively. We assumed here 
that l3Np e N. 

As £ n- V'^{e) is concave in Regime b for all i e N, if the elongation of a 
spring in the local minimum is altered by an arbitrarily small perturbation, 
it will move (according to the sign of the perturbation) to either or . As 
a consequence, any system of N springs with q is unstable and we may in 
the following restrict to the case q = 0. 

The average strain of a system with l3N springs in equilibrium and the first 
PNp springs in phase a fulfills the identity 

^ I3N ^ PN 1 ^ 



av y i ^ ^ * ' TV ^ S + a' 

i=l i=l3Np+l i=l3N+l 

Here we only study the limiting case (5 \, where there is no further small 
correction of the damaged springs. In this case, we obtain 

KQ+a('m) Kq + a(m) 

where we used (Al)-(A4), ^ and ([S]); especially ej. = + e*. 

Resolving ([7]) , we get the stress-strain relationship for a TV-springs system 

ko{e - Ep + P - 1) 

with Sp{m) := — p)£*. The latter can in a natural way be identified with 
the plastic strain. From we see that a only depends on m and on the elastic 
strain Sd :—£ — £„. 



2.3 Analytic computation of the continuum limit 

We identify the continuum limit of in the framework of F-convergence (see, 
e-g-i [4j). This is not a standard procedure as V also depends on the spatial 
position. 

Let L = 1 (which can always be obtained by rescaling), and := (0, 1). For 
prescribed / > 0, we impose the boundary conditions 

u^ = 0,u% = l. (9) 

This selection of boundary conditions appears natural for the discrete system 
in the absence of next-to-nearest neighbor-interactions. 

We now specify our assumptions on V. The mechanical pair-potential is a 
function of the deformation gradient. In addition, in order to capture the effect 
of damage expressed in the assumptions (A1)-(A4) above, we need to respect 
the dependence of V on the spatial position i. Let dom(F) = [0, 1] x D for 
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suitable D C K be the domain of definition of V. We postulate the following 
conditions on V. 

(COND 1) There exist positive constants ci, Ci, such that 

ci(|e| - 1) < Vix,e) < Ci{\e\ + 1) for all {x,e) e [0,1] x D. 

(COND 2) The function x n- V{x,e) is continuous for any e E D. 

After extending V given by ^ continuously in i by interpolation, we easily verify 
that this extension (again called V) satisfies both COND 1 and COND 2. For 
given deformations un ■— {uf^)o<i<N, we may rewrite the overall mechanical 
energy of the chain as a functional of the discrete displacement gradient, 

JV-l . N _ N 

E^M J2 ^(^' "'^fe ) • (10) 

2 = ^ 

The functional V coincides with ([2]), but the dependence on the subscript i has 
moved to the first argument. We scaled En by /iat as we are dealing with a 
microscopic energy. 

Following ideas in [5], we introduce for N <E N the set of all functions 
u : /iatZ n [0, 1] R, setting Ui :— u{ihN)- We tacitly identify -4(0, 1) with the 
piecewise affine linear interpolations, i.e. 

A{id, 1) {m : [0, 1] -> R I w is affine in {ihN, (i + l)/iAr), Q<i<N-l]. 

With Ui = u{ihN), as a shorthand notation for the second argument in ()10|) . we 
introduce the symbol 

V^un{x) X e [ihN, {i + l)hN), Q<i<N -I, 

hN 

which is a discrete approximation of Vu with step size /ijv. 
For later use we introduce E'^ : ^(0, 1) — R U {+00} by 



ENiu), if m(0) = 0, u(l) = /, 
+00, otherwise. 



Theorem 1 ('Continuum limit of En) 

Let V satisfy the conditions ( COND 1), ( COND 2) and let I > 0. Then the func- 
tional E^j^ converges in the T -sense for N — > 00 a functional E^ : i^(0, 1) — >■ M 
defined by 

EUuV.J ^o*ix, V7/(a:)) dx ifue i7i^i(0, 1), u(0) = 0, u{l) = /, 
\ +c>o else. 

Here, Vo(a;, z) :— i min |p^(a;, zi) + V{x,Z2) \ zi + Z2 = 2z|, and Vq* denotes 
the convexification of Vo (see, e.g., [17 ). 

Proof (a) Proof of the lim inf -inequality. 
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Let a sequence (wAr)7veN C L^(0,1) be given with ujv u in L^{0,1) as 
iV — T' oo. W.l.o.g. E'-j^{un) < oo and thus En{un) < oo. Now, as the bounds 
in (COND 1) are uniform in x, we can apply Theorem 2 in [2] which proofs the 
lim inf-inequality for (ujv)Ar- 

(b) Proof of the limsup-inequality. 

Let u € L^(0,1) be given. We have to show the existence of a sequence 
{un)n£N C L^{0, 1) such that un ^ u in L^{0, 1) as TV ^ cx) and 



By a standard density argument, cf. Part (b3) in the proof of Theorem 2 in [5], 
we may restrict to the afHne case u{x) = ax + b. 

For the construction we first ignore the boundary conditions (O. Let N = 
km for some m, k (z N. We will choose functions un which are periodic in any 
subinterval of (0, 1) with length mh^- We write the integrand as V{dM{x),-)^ 
where dAr:(0,l)^Kisa piecewise constant bounded function with djv — >■ Id 
for N ^ oo. The method can be generalized to more general situations. With 
these settings we find 



where (f := dN{{i + 1/2)to/ijv); and (COND 2) has been used. 

By convexity of Vq* and Caratheodory's theorem (see, e.g., [T7]) we know 
that for any < i < k — 1, there exist real numbers Ai, A2 with < Ai, A2 < 1 
and < Ai + A2 < 1 such that 

V„**{d\a) - XiVo{d^'\a') + X2Vo{d^'\a^) + {1 - Xi - X2)Void^^a^), (12) 
{d\a) = id\Vu) = Al(rf''^al) + A2(d*'^a2) + (l-Al-A2)((i*'^a3). (13) 

For given Ai, A2 we introduce the sets 

rjf := f^nu(^Z"o {ihm,ihm + h[Xim\], 

:= 17 n Ufjo {{ihm + h[Xim\), (ihm + h[{Xi + A2)toJ], 
-.^nn U^-^ {{ihm + h[{Xi + A2)mJ), {i + l)hm], 

such that fl^nf' Ufl^ un^. 

By definition of Vq it holds for < i < fc — 1 and any s = 1, 2, 3, 



limsupi?]v(uAf) < E {u). 




(11) 
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We choose un{x) ~ a]^{x)x + &, where 



ajv(a;)=< 
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{{im+2j)hN, ini + 2j + l)hN], 
{{iin + 2j)hN, {ini + 2j + l)h]y], 
{{im+2j)hN, {im + 2j + l)h]y], 
{{im+2j + l)hN, iim+2j + 2)hN], 
{{im + 2j + l)hN, iim+2j + 2)hN], 
{{im + 2j + l)hN, {im+2j + 2)hN]. 



For s = 1,2,3 and < i < fc — 1 we set 

' df, if a; e r^f n UfJ^^^ {{im + 2j)hN, {im + 2j + l)h 



iixen^ n □"f'o ^ {{im + 2j + l)hN, (im + 2j + 2)h 



The ansatz for depends on the current interval, as V is x-dependent. 
Eqn. now reads after setting [A3TOJ :—m— [AimJ — [A2mJ for short 



k-i 3 



1 r. 



EM{uN) = hNY.Y.\-^^'^^ 9 V{d'{\al) + V{d 



2 I "2; 



1=0 s=l 
fe-1 3 



(14) 



^^^^^^ 

i=Q s=l 

For m — >■ 00 we have [A^mJ/m — As, s = 1, 2. Consequently, using (jl2p . 

En{un) ^ [ V„**{x,\7u{x))dx = E\u) as TV ^ 00. 

We still have to show that ujy ^ u in £^(0, 1). If in the derivation of Eqn. (|14l) 
we formally set l^(a;, v) := v (which is feasible), we obtain 

r-l 



V.^(.)dx = i^ai + i^a^ 



m 

00 as above 



1 



[AitoJ [A2TOJ 



and in the limit m 

lim / VM7v(a^) dx = Aia^ + A2a^ + (1 — Ai — A2)a^ 



Vu{x)dx, (15) 



where Eqn. (fT3| has been used. By (|15p . mat ^ u in L''(0, 1) for 1 < r < cxd. 

We still need to incorporate the boundary condition ([9]) . If mat is the recovery 
sequence from above, we define 

r ^^^x, xe[o,h^), 

vn{x) := < un{x), X <E [hN,l - ^n], 

[ '~^-l^l"''-\ x~l) + l, xe{l-hr,,l]. 

By construction, vn is continuous on [0,1] with vn{0) — 0, WAr(l) = I, and 
lim ff-^ 00 {En{un) — En{vn)) — 0, which shows that ujv is the sought recovery 
sequence. □ 
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3 Numerical results 



We study in this section the numerical convergence of the stress-strain response 
of finite mass-spring systems to the continuum hmit of Eqn. ([8]). We analyze the 
overall loading-unloading response of discrete systems composed oi N = x + y 
springs, where x denotes the number of undamaged springs, while y specifies the 
number of damaged springs. For the spring constants, we use the parameters 
ko = 50.00 X 10^ Pa, h = -22.44 x 10^ Pa, 0.25, e., = 0.52, 6 = 50.00 x 10^ 
Pa, which correspond to aa = 16.67 MPa, Cc = 5.00 MPa, and Act = CTc ^ fa = 
-11.67 MPa. 



0.5r 0.5r 




Figure 1: Numerical overall stress-strain response of bistable mass-spring chains 
for different total number of springs A'^. 



Fig. [T] shows the results for different x, y. The stress-strain curves of the 
discrete microscopic chains follow sawtooth patterns which alternate elastic and 
plastic steps. The plastic steps are characterized by fixed total strain and micro- 
scopic branch switching. Such sawtooth-like responses converge to a perfectly 
plastic behavior with a loading (cr = (Xa) and an unloading plateau (u = cTc) for 
increasing values of N, as predicted by Eqn. ([8]). It is worth observing that the 
stress is zero for e < 1 — /3, while for e > 1 — /3 that the system is able to bear 
stresses tr > 0. Fig. [T] shows, in addition, that the energy dissipation capacity 
of the system (area enclosed by stress-strain curve in a loading-unloading cy- 
cle) reduces for decreasing /?, that is for increasing number of damaged springs 
(bottom-right panel). Remarkably, the quantity 1-/3 can be regarded as the 
'activation' strain of the system. 
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4 Concluding remarks 



In this paper, for the first time a bistable-spring ansatz has been proposed 
capable of incorporating damage in plastically deformed materials. 

We have shown that a suitable modification of the model recently proposed 
by [S] for CNT foams is able to handle preconditioning induced material damage, 
which is characterized by an activation strain different from zero; a reduction 
in the energy dissipation capacity; and permanent deformation. The latter, in 
particular, coincides with the activation strain. The new model allows us to 
extend the 'transformational plasticity' concept discussed by |15j from time- 
independent hysteretic behavior to fatigue-type material damage. It applies to 
a wide class of materials showing Mullins-like behavior, [13], which includes 
besides CNT arrays rubber-like and soft biological materials. 

We address a multiscale formulation of the present model, accounting for 
graded mechanical properties along the height of the structure, [16], and the 
analytic computation of the energy dissipated by the system in the continuum 
limit in future work. 
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